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Abstract 

Background: Rabies virus (RABV) causes a fatal infection of the central nervous systems (CNS) of warm-blooded 
animals. Once the clinical symptoms develop, rabies is almost invariably fatal. The mechanism of RABV pathogenesis 
remains poorly understood. Recent studies have shown that microRNA (miRNA) plays an important role in the 
pathogenesis of viral infections. Our recent findings have revealed that infection with laboratory-fixed rabies virus 
strain can induce modulation of the microRNA profile of mouse brains. However, no previous report has evaluated 
the miRNA expression profile of mouse brains infected with RABV street strain. 

Results: The results of microarray analysis show that miRNA expression becomes modulated in the brains of mice 
infected with street RABV. Quantitative real-time PCR assay of the differentially expressed miRNAs confirmed the 
results of microarray assay. Functional analysis showed the differentially expressed miRNAs to be involved in many 
immune-related signaling pathways, such as the Jak-STAT signaling pathway, the MAPK signaling pathway, 
cytokine-cytokine receptor interactions, and Fc gamma R-mediated phagocytosis. The predicted expression levels of 
the target genes of these modulated miRNAs were found to be correlated with gene expression as measured by 
DNA microarray and qRT-PCR. 

Conclusion: RABV causes significant changes in the miRNA expression profiles of infected mouse brains. Predicted 
target genes of the differentially expression miRNAs are associated with host immune response, which may provide 
important information for investigation of RABV pathogenesis and therapeutic method. 

Keywords: Street strain rabies virus, Brain infection, MicroRNA profiling, Gene profiling, Target prediction, Functional 
enrichment 



Background 

The rabies virus (RABV), a member of the family Rhab- 
doviridae, is a highly neurotropic virus that can cause 
fatal infections of the central nervous systems (CNS) of 
warm-blooded animals [1,2]. Although significant advances 
have been made in rabies prevention and control, the 
disease remains a major threat to public health. It 
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causes 55,000 people die around the world every year [3]. 
Despite the catastrophic clinical outcome of RABV en- 
cephalomyelitis, the histopathological changes observed in 
the CNS are typically relatively mild, showing varying 
degrees of mononuclear inflammatory cell infiltration of 
the leptomenings, perivascular cuffing, microglial activa- 
tion, and neuronophagia. Although there are several hy- 
potheses under active study at present, the pathogenesis of 
the rabies virus has not yet been determined. 

MiRNA (miRNA, miR) is endogenous w 22 nt RNA 
that negatively regulates gene expression by translational 
repression [4], It binds to the complementary sequences 
in the mRNAs and blocks the translation or accelerates 
mRNA decay [5]. MiRNAs play key roles in cellular 
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processes such as development, differentiation, cell pro- 
liferation, and hematopoiesis [6-9]. Recently, evidence 
has demonstrated that cellular miRNAs exert regulatory 
functions in virus-host interactions [10,11]. It is becom- 
ing increasingly clear that miRNAs of cellular origin can 
positively or negatively influence viral infection. For ex- 
ample, miR-122 is indispensable to replication of the 
hepatitis C virus (HCV), whereas miR-196 and miR-296 
substantially attenuate viral replication [12,13]. A recent 
study reported that miR-28, miR-125b, miR-150, miR- 
223, and miR-382 inhibit replication of the human im- 
munodeficiency virus (HIV) in CD4 + T cells [14]. 

Microarray analyses have been recently employed to 
detect changes in host miRNA expression, which can help 
reveal molecular pathways that govern viral pathogenesis. 
By using miRNA microarray profiling, researchers have 



observed differentially expressed patterns of cellular miR- 
NAs in the lungs of mice infected with influenza virus 
[15]. Another study found miRNAs to be significantly 
regulated in mouse brains upon Venezuelan equine en- 
cephalitis virus (VEEV) infection [16]. Our own recent 
findings suggest that infection with laboratory-fixed rabies 
virus, ERA (Evelyn Rokitnicki Abelseth), can induce 
modulation of the microRNA profile of the mouse brain 
[17]. However, no report has yet been made regarding the 
assessment of the host miRNA expression profile in 
mouse brains upon infection with the street strain of 
RABV. 

In this study, we performed an expression profile of 
cellular miRNAs in the brains of mice infected with the 
highly pathogenic street rabies virus. Meawhile,we per- 
formed target prediction and functional enrichment of 
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Figure 1 Outcomes of mice infected with RABV Fujian strain. After i.e. injection of 10 5 ffu of street Fujian RABV strain, (A) clinical score and 
(B) copy number of RABV N mRNA were recorded as described in Materials and Methods. Mice were monitored for survival for 21 days. Data 
were obtained from 8 mice (three mice for N mRNA) in each group. Data are the mean ± standard deviation (SD) of one representative 
experiment. Similar results were obtained in three independent experiments. 
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the differentially expressed miRNAs. It was shown that 
several miRNAs were modulated in mouse brains infected 
with RABV. Finally, we performed gene microarray ana- 
lysis and qRT-PCR measurement to verify the expression 
levels of the predicted targets of the modulated miRNAs 
in these pathways. The results of functional enrichment 
revealed that many of the predicted targets of these miR- 
NAs play key roles in the immune response, which are 
known to be associated with the pathogenesis of RBAV. 



Results 

Characterization of pathogenicity of RABV Fujian strain 
in mice 

All infected animals showed RABV-specific symptoms 
that increased in severity in a time-dependent manner. 
After inoculation, all mice developed the clinical signs of 
disordered movement at 4 days post-infection (dpi). At 6 
dpi, considerable aggravation of typical clinical signs 
was observed, with the onset of trembling, shaking, 
anger, and hyperexcitation followed by general paralysis 
(Figure 1A). As presented in a previous study, out of all 
infected animals, 12.5%, 25%, and 50% were dead at 6, 7, 
and 8 dpi, respectively, and all mice succumbed to RABV 
at or before 9 dpi [18]. Viral load in the brain was 



monitored up to 9 dpi using Taqman qRT-PCR 
(Figure IB). RABV replicated rapidly in brain with an in- 
crease in copy number from 3 dpi and reached a maximal 
viral load at 7 dpi. The results demonstrate that the RABV 
Fujian strain is highly pathogenic in mice. 



Modulation of miRNA profile in brain in response to 
RABV infection 

To determine changes in miRNA expression in mouse 
brains in response to street RABV infection, we evaluated 
miRNA expression profiles at 7 dpi. The two-way hierarch- 
ical cluster heat map clear showed different expression pat- 
tern of host miRNAs between RABV and mock infections 
(Figure 2A). MiRNAs whose relative expression levels 
showed a fold change (FC) > 2 and P < 0.01 were consid- 
ered significantly up-regulated, and those with FC < -2 and 
P<0.01 were considered significantly down-regulated. As 
shown in Figure 2B, nine miRNAs, miR-691, miR-377, 
miR-1935, miR-190, miR-1902, miR-135a*, miR-203, miR- 
2138, and miR-290-5p, were found to be significantly up- 
regulated. However, only one miRNAs, miR-145, was 
found to be down-regulated upon RABV infection. This 
indicates that host miRNAs were modulated in the CNS 
upon infection with street rabies virus. 
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Figure 2 MiRNA profile of street RABV-infected mouse brain. (A) Two-way hierarchical cluster heat map showing all significantly expressed 
miRNAs in three independent samples (P<0.01). Each row shows the relative expression level of a single miRNA. Each column shows the 
expression level of a single sample. Up-regulated miRNAs are shown in red and down-regulated miRNAs are shown in green. Significantly 
differentially expressed miRNAs in mouse brain upon RABV infection by microarray analysis. MiRNAs whose relative expression levels showed a 
fold change (FC)>2 and P < 0.01 were considered significantly up-regulated, and those with FC<-2 and P < 0.01 were considered significantly 
down-regulated. 
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Confirmation of differentially expressed miRNAs by 
qRT-PCR 

To validate the differential expression profiles of miRNAs 
obtained by microarray analysis, quantitative RT-PCR was 
performed on six selected differentially expressed miRNAs 
including miR-691, miR-377, miR-1935, miR-190, miR- 
203, and miR-145. The data demonstrate that the overall 
results of qRT-PCR were consistent with those of the 
microarray analysis. Although differences were observed 
between these two types of analysis due to intrinsic dif- 
ferences between the techniques, the qRT-PCR results 
showed the same relative regulation of differentially 
expressed miRNAs as the microarray data results 
(Figure 3). 

Target prediction and functional analysis of differentially 
expressed miRNAs 

Target genes regulated by these differentially expressed 
miRNAs were predicted using Targets can Mouse, 
MicroCosm, and miRanda. For these ten differentially 
expressed miRNAs, TargetScan predicted 2,058, Micro- 
Cosm predicted 5,433 and miRanda predicted 29,742 
target genes. Of these, 3,038 target genes were predicted 
under all three systems (Additional file 1: Figure SI). 
Gene ontology (GO) analysis in the Database for Anno- 
tation, Visualization and Integrated Discovery (DAVID) 
was performed for these miRNAs using the predicted 
gene targets [19]. Functional analysis revealed 106 GO 
terms to be involved in biological processes, 14 in 




microRNAs 

Figure 3 Verification of differentially expressed miRNAs by 
qRT-PCR. Six differentially expressed miRNAs were selected from 
miRNA microarray datasets and examined by qRT-PCR. The fold 
change from the qRT-PCR was determined using the 2 _AACt method 
and all miRNA expression values were normalized against the U6 
endogenous control. Data from qRT-PCR are shown as 
mean ± standard deviation (SD) of one representative experiment. 
Similar results were obtained in three independent experiments. 
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molecular function, and 20 in cellular components 
(P < 0.01) (Additional file 2: Table SI). The twenty most 
common GO categories were cellular processes, meta- 
bolic processes, cellular metabolic processes, macromol- 
ecule metabolic processes, and cellular macromolecule 
metabolic process (Figure 4). These analyses suggest that 
cellular miRNAs may regulate cellular metabolic pro- 
cesses during street RABV infection, either directly or 
indirectly. 

Pathway analysis of target genes of differentially 
expressed miRNAs 

To identify the biological pathways that become active 
in the mouse brain in response to RABV infection, we 
mapped the target genes of differentially expressed miR- 
NAs to canonical signaling pathways in the Kyoto 
Encyclopedia of Genes and Genomes (KEGG). The 
results showed that 25 statistically remarkable categories 
(P < 0.05) were enriched (Additional file 3: Table S2). As 
shown in Table 1, the predicted target genes of six up- 
regulated miRNAs, miR-691, miR-377, miR-1935, miR- 
190, miR-203, and miR-135a*, and one down- regulated 
miRNA, miR-145, were found to be involved in 
immune-related pathways, such as the Jak-STAT signal- 
ing pathway, MAPK signaling pathway, Fc gamma R- 
mediated phagocytosis and cytokine- cytokine receptor 
interactions. The predicted target genes of five up- 
regulated miRNAs, miR-691, miR-377, miR-190, miR-203, 
and miR-1290-5p, and one down-regulated miRNA, miR- 
145, were found to be involved in other pathways, such as 
the Adherens junction, Wnt signaling pathway, Axon 
guidance, cell cycle, TGF-beta signaling pathway, and 
Focal adhesion. 

DNA microarray assay and qRT-PCR measurement of 
miRNA targets 

MiRNAs predominately function as repressors of target 
gene expression. The miRNAs and their targets show 
mutually antagonistic expression levels. To determine 
whether any such correlation exists between deregulated 
miRNA levels and their corresponding targets, we per- 
formed DNA microarray assay and qRT-PCR validation. 
To identify the genes involved in the pathways common 
to miRNA target prediction and RABV infection, we 
measured the expression of genes from Jak-STAT signal- 
ing pathway (SOCS4), cytokine-cytokine receptor inter- 
actions (7125, CD40, VEGFA, and CCR10), MAPK 
signaling pathway (MAPKAPK3), Fc gamma R-mediated 
phagocytosis (ARF6), and Wnt signaling pathway 
(NFAT5). As presented in Table 2, MAPKAPK3 and 
IL25, the targets of up-regulated miRNA miR-691; 
SOCS4, the target of up-regulated miRNAs miR-377; 
CCR10 and NFAT5, the targets of up-regulated miRNA 
miR-1935; VEGFA, the target of up-regulated miRNA miR- 
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GO:0010468 regulation of gene expression 
GO:0031326 regulation of cellular biosynthetic process 
GO:0009889 regulation of biosynthetic process 
GO:0007275 multicellular organismal development 
GO:0080090 regulation of primary metabolic process 
GO:0060255 regulation of macromolecule metabolic process 
GO:0006139 nucleobase, nucleoside, nucleotide and nucleic acid metabolic. 

GO:0044249 cellular biosynthetic process 
GO:0031323 regulation of cellular metabolic process 
GO:0032502 developmental process 
GO:0009058 biosynthetic process 
GO:0019222 regulation of metabolic process 
GO:0034641 cellular nitrogen compound metabolic process 
GO:0006807 nitrogen compound metabolic process 
GO:0044260 cellular macromolecule metabolic process 
GO:0043170 macromolecule metabolic process 
GO:0044237 cellular metabolic process 
GO:0044238 primary metabolic process 
GO:0008152 metabolic process 
GO:0009987 cellular process 
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Figure 4 Enriched GO terms in the biological process category among differentially expressed miRNAs. After miRNA microarray 
assay, significantly enriched GO analysis in the biological process category was performed on differentially expressed genes in the brains of 
RABV-infected mice using DAVID (P< 0.01). Only the top twenty GO terms are listed here. For other enriched GO terms, please see 
Additional file 1: Table SI. 
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203; and CD40, the target of up-regulated miRNA miR- 
290-5p, were found to be down-regulated. ARF6, the 
target of down-regulated miRNA miR-145, was found 
to be up-regulated. The negative correlation between 
these miRNAs and their targets showed these pathways 
to be part of the RABV infection response at the expres- 
sion level Together, these results strongly suggest that 
certain miRNAs may be associated with RABV infection 
and pathogenesis. 

Discussion 

RABV, a pathogen well-adapted to the nervous system, 
infects the neurons of warm-blooded animals. Despite 
the catastrophic clinical outcome of RABV encephalo- 
myelitis, especially that caused by street viruses, the 
histopathological changes observed in the CNS are typic- 
ally relatively mild. CNS has more intrinsic mechanisms 
for controlling immune response than other organs. In- 
deed, many viral infections can be cleared from the CNS 
by immune mechanisms [20]. This may indicate that, 
upon viral infection, the immune privilege of the CNS is 
not as strong as has been proposed [21]. 

In mammals, approximately 30% of all protein-coding 
genes are predicted to be regulated by miRNAs [22]. 
Currently, nearly a thousand miRNAs have been cloned, 
each potentially regulating hundreds of genes by 
complementary binding to the 3 '-untranslated region 
(3'-UTR) of the target mRNAs. Recent publications have 



provided compelling evidence that a range of miRNAs 
are involved in the regulation of immunity, including the 
development and differentiation of B and T cells, prolif- 
eration of monocytes and neutrophils, antibody switch- 
ing and the release of inflammatory mediators [4,23-25]. 
More importantly, researchers have reported that cellu- 
lar miRNAs play key regulatory roles during viral infec- 
tion and that altered cellular miRNA expression in 
response to viral infection may be an important deter- 
minant of virulence [10,11,26]. 

Recently, our findings suggested that infection with 
laboratory-fixed rabies virus strain, ERA, induced modu- 
lation of microRNA profile of mouse brains [17]. In the 
current study, a comprehensive examination of miRNA 
expression from brains of street-RABV- and mock- 
infected mice was performed. The results showed that 
host miRNAs were modulated upon RABV infection. 
We found that the expression profiles of host miRNAs 
in mouse brains infected by these two strains of RABVs 
were completely different from each other. Interestingly, 
the same differences in transcriptome were observed be- 
tween the mouse brains upon infection with the two 
strains of RABVs, ERA and street RABV Fujian strain 
(unpublished data) [18]. This discrepancy can be partly 
attributed to differences in virulence of the viruses. Fur- 
ther study is required to determine the correlation be- 
tween specific miRNA expression and RABV virulence 
using reverse genetics or RNA interference. Some 
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Table 1 Predicted targets of modulated miRNAs upon RABV infection involved in immune response pathways 



miRNA 


KEGG pathway 


Target genes 


Number of 
targets 


P value 


mmu-miR-377 


mmu04630: Jak-STAT signaling pathway 


CSF2RB2, EP300, JAK2, SOCS4, BCL2L 


5 


0.01 


mmu-miR-377 


mmu04010: MAPK signaling pathway 


CACNA2D1, FASL, PPM] A, PPP3R 1, RASA 1, FLNC, PLA2G4A, SRF 


8 


0.02 


mmu-miR-377 


mmu04060: Cytokine-cytokine receptor 
interaction 


CSF2RB2, FASL, ACVR2A, VEGFA, IL18RAP, KITL, PDGFRA, LIB 


8 


0.001 


mmu-miR-377 


mmu04666: Fc gamma R-mediated 
phagocytosis 


WASF1, LIMK1, PRKCD, PLA2G4A, RPS6KB1 


5 


0.01 


mmu-miR-377 


mmu04520: Adherens junction 


M//\CX7 PDinn VCC 7 D\/DI A DTDDQ QAAAH/1 
VVAjrl, trjUU, Ytjl, rvKLH, rlrnb, jIVIAUh 


0 


I . I ot-Uo 


mmu-miR-377 


mmu04310: Wnt signaling pathway 


PDinn \A/KITQA DDD3D1 C7D3 P7D/I AYIM1 RTDr CA/MD/f 
trjUU, VvNljA, rrrjnl, rZUj, r/LUH, AAiN / , D/nL, J/VIAU4-, 

CUL1, DKK1 


1 n 
I U 


o.Dzt-U4 


mmu-miR-377 


mmu04360: Axon guidance 


DDD1D1 DAQA1 1 IhAI/l QCAAA^A C/TA/M3P) 

rrrjnl, HAjAI, LIIVIM, jcIVIAoA, jLIVIAjU 


r 
J 


U.UU I 


mmu-miR-377 


mmu041 10: Cell cycle 


EP300, ESPL1, CDC26, CUL1 


4 


0.007 


mmu-miR-377 


mmu04350: TGF-beta signaling pathway 


EP300, ACVR2A, THBS 1, PITX2, SMAD4, CULT, LEFTY2, RPS6KB1 


8 


6.96E-06 


mmu-miR-377 


mmu04510: Focal adhesion 


MYLK2, XIAP, ITGA6, VEGFA, THBS1, PDGFRA, FLNC, LAMC1 


8 


0.018 


mmu-miR-691 


mmu04630: Jak-STAT signaling pathway 


GHR, SPRY2, TYK2, SPRED2, IL12A, CLCF1 


6 


0.01 


mmu-miR-691 


mmu04010: MAPK signaling pathway 


PLA2G10, MKNK2, MAPKAPK3, CACNB2, TGFBR2, DUSP9 


5 


0.02 


mmu-miR-691 


mmu04060: Cytokine-cytokine receptor 
interaction 


GHR, ACVR2A, BMP2, IL25, AMHR2, IL12A, TGFBR2, TNFSF14, CLCF1 


9 


0.001 


mmu-miR-691 


mmu04666:Fc gamma R-mediated 
phagocytosis 


IVM /A/DDCPl 
LYN, INrrDU 


z 


U.U I 


mmu-miR-691 


mmu04520: Adherens junction 


CSNK2B, SRC, CDH1, TGFBR2 


4 


1.18E-06 


mmu-miR-691 


mmu04310: Wnt signaling pathway 


DKK2, CSNK2B, CAMK2G, CTBP2, VANGL1 


5 


3.62E-04 


mmu-miR-691 


mmu04360: Axon guidance 


DTIS~> CEAAAtZD AQIIAA1 MDD1 Dl VKIA 1 CD/~/\D3 QEhAAZA 

rll\z, jLIVIAoU, AdLIIVII, Nnrl, rLANAz, bnuArj, jcIVIAdA 


7 


0.001 


mmu-miR-691 


mmu041 10: Cell cycle 


YWHAB, TFDP1 


2 


0.007 


mmu-miR-691 


mmu04350: TGF-beta signaling pathway 


ACVR2A, BMP2, THBS1, AMHR2, TFDP1, SP1, TGFBR2 


7 


6.96E-06 


mmu-miR-691 


mmu04510: Focal adhesion 


PTK2, IBSP, ITGA10, PGF, CAPN2, THBS 1, SRC 


7 


0.018 


mmu-miR-1935 


mmu04630: Jak-STAT signaling pathway 


LIFR, SPRY3 


2 


0.01 


mmu-miR-1935 


mmu04010: MAPK signaling pathway 


MAP3K12, FLNC, MART 


3 


0.02 


mmu-miR-1935 


mmu04060: Cytokine-cytokine receptor 
interaction 


ACVR1, HGF, LIFR, CCR10 


4 


0.001 


mmu-miR-1 935 


mmu04666: Fc gamma R-mediated 
phagocytosis 


l\AAQCI/Q 
IVIAr\L.I\j> 


i 


n m 

U.U I 


mmu-miR-1935 


mmu04520: Adherens junction 


YES1, PVRL1 


2 


1.18E-06 


mmu-miR-1 935 


mmu04310: Wnt signaling pathway 


NrA 1 D 


1 


3.62E-04 


mmu-miR-1935 


mmu04360: Axon guidance 


NFAT5, PLXNA4, SRGAP2, DCC 


4 


0.001 


mmu-miR-1935 


mmu041 10: Cell cycle 


STAG1, CDKN1C 


2 


0.007 


mmu-miR-1935 


mmu04510: Focal adhesion 


FLNC 


1 


0.018 


mmu-miR-1 90 


mmu04630: Jak-STAT signaling pathway 


IL6 


1 


0.01 


mmu-miR-1 90 


mmu04010: MAPK signaling pathway 


FASL, PLA2G6, BDNF, FGF14 


4 


0.02 


mmu-miR-1 90 


mmu04060: Cytokine-cytokine receptor 
interaction 


FASL, HGF, IL6, CXCL16 


4 


0.001 


mmu-miR-1 90 


mmu04666: Fc gamma R-mediated 
phagocytosis 


PLA2G6, ARPC5 


2 


0.01 


mmu-miR-1 90 


mmu04520: Adherens junction 


SMAD2, TJP1 


2 


1.18E-06 


mmu-miR-1 90 


mmu04310: Wnt signaling pathway 


SMAD2, PPP2R1B, ROCK1, PLCB1 


4 


3.62E-04 


mmu-miR-1 90 


mmu04360: Axon guidance 


EPHA1, ROCK! 


2 


0.001 


mmu-miR-1 90 


mmu041 10: Cell cycle 


SMAD2, CHEK1, RBL1 


3 


0.007 
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Table 1 Predicted targets of modulated miRNAs upon RABV infection involved in immune response pathways 

(Continued) 



mmu-miR-190 


mmu04350: TGF-beta signaling pathway 


SMAD2, PPP2R1B, PITX2, ROCK1, RBL1 


5 


6.96E-06 


mmu-miR-190 


mmu04510: Focal adhesion 


CAPN2, HGF, R0CK1 


3 


0.018 


mmu-miR-135a* 


mmu04630: Jak-STAT signaling pathway 


EP300, IL24, IFNAR2, IL 12RB 1JFNG 


5 


0.01 


mmu-miR-135a* 


mmu04060: Cytokine-cytokine receptor 
interaction 


IL24JFNAR2JL 12RB 1,IFNG,CXCL 16,CD40LG 


6 


0.001 


mmu-miR-135a* 


mmu04666: Fc gamma R-mediated 
phagocytosis 


NCF1 


1 


0.01 


mmu-miR-135a* 


mmu04520: Adherens junction 


EP300, ACTN2 


2 


1.18E-06 


mmu-miR-135a* 


mmu04310: Wnt signaling pathway 


EP300, WNT1, PPP2R1A 


3 


3.62E-04 


mmu-miR-135a* 


mmu041 10: Cell cycle 


EP300, YWHAQ, ANAPC7 


3 


0.007 


mmu-miR-135a* 


mmu04350: TGF-beta signaling pathway 


EP300, BMP4, IFNG, PPP2R1A 


4 


6.96E-06 


mmu-miR-135a* 


mmu04510: Focal adhesion 


PTEN, ITGA8, ACTN2 


3 


0.018 


mmu-miR-290-5p 


mmu04630: Jak-STAT signaling pathway 


IL13RA1, IFNGR2, AKT1, CSF2RB 


4 


0.01 


mmu-miR-290-5p 


mmu04010: MAPK signaling pathway 


GADD45A, MAPKAPK3, MEF2C, MAP3K1, AKJ1FGF23 


5 


0.02 


mmu-miR-290-5p 


mmu04060: Cytokine-cytokine receptor 
interaction 


IL13RA1, IL18R1, CD40, IFNGR2, CXCL12, ACVR1B, CSF2RB 


7 


0.001 


mmu-miR-290-5p 


mmu04666: Fc gamma R-mediated 
phagocytosis 


MARCKS, AKJ1 


2 


0.01 


mmu-miR-290-5p 


mmu04520: Adherens junction 


LM07, ERBB2, PTPRB, IFF), INSR, SSX21P, PVRL3 


7 


1.18E-06 


mmu-miR-290-5p 


mmu04310: Wnt signaling pathway 


IFF], CAMK2G, PLCB1, AXIN2 


4 


3.62E-04 


mmu-miR-290-5p 


mmu04360: Axon guidance 


ROBOT, L1CAM, SEMA5B, CXCL12 


4 


0.001 


mmu-miR-290-5p 


mmu041 10: Cell cycle 


SMC3, MCM5, GADD45A, CDC25A, TFDP1, HDAC2, RAD21, ORC5L 


8 


0.007 


mmu-miR-290-5p 


mmu04350: TGF-beta signaling pathway 


TFDP1 


1 


6.96E-06 


mmu-miR-290-5p 


mmu04510: Focal adhesion 


XIAP, ERBB2, AKF1 


3 


0.018 


mmu-miR-203 


mmu04630: Jak-STAT signaling pathway 


1124, CNTFR, IL22RA2, CCND1, S0CS3, AKJ2, IL12B 


7 


0.01 


mmu-miR-203 


mmu04010: MAPK signaling pathway 


CACNA2D1, MAP3K13, DUSP5, MAP4K3, STK3, NLK, AKT2, ATF2, 
PDGFRA, RAP1A, MAPI MAP3K1, PPM1B, FGF7, CRK, PRKCA 


16 


0.02 


mmu-miR-203 


mmu04060: Cytokine-cytokine receptor 
interaction 


MET, 1124, CNTFR, ACVR2A, IL22RA2, VEGFA, XCL1, PDGFRA, IL12B 


9 


0.001 


mmu-miR-203 


mmu04666: Fc gamma R-mediated 
phagocytosis 


PLD2, AKF2, DNM1L, PPAP2B, CRK, VAV3, PRKCA 


7 


0.01 


mmu-miR-203 


mmu04520: Adherens junction 


MET, NLK, SRC, PTPRB, SNAI2 


5 


1.18E-06 


mmu-miR-203 


mmu04310: Wnt signaling pathway 


SFRP2, CCND1, NLK, APC, CSNK1A1, CUL1, PRKCA 


7 


3.62E-04 


mmu-miR-203 


mmu04360: Axon guidance 


MET, ABL1, SEMA3A, SEMA6A, SEMA5A 


5 


0.001 


mmu-miR-203 


mmu041 10: Cell cycle 


YWHAQ, CCND1, ABL1, SKP2, CUL1 


5 


0.007 


mmu-miR-203 


mmu04350: TGF-beta signaling pathway 


ACVR2A, ID2, BMP5, CUL1, ID4 


5 


6.96E-06 


mmu-miR-203 


mmu04510: Focal adhesion 


MET, PXN, CCND1, VEGFA, PPP1C, SRC, PDGFRA, RAP1A, RAPGEF1, 
TNC, COL4A4, CAV1, CRK, VAV3, PRKCA 


15 


0.018 


mmu-miR-145 


mmu04630: Jak-STAT signaling pathway 


MYC, STAM, STAT4, S0CS7 


4 


0.01 


mmu-miR-145 


mmu04010: MAPK signaling pathway 


MAP3K3, HSPA1L, PPP3CA, DUSP4, MYC, RASA2, RAPGEF2, FLNB, 
DUSP6, TGFBR2, MAP4K4, CRKL 


12 


0.02 


mmu-miR-145 


mmu04060: Cytokine-cytokine receptor 
interaction 


FLT3L, TNFRSF1 IB, IL17RB, INHBB, ACVR1B, TGFBR2 


6 


0.001 


mmu-miR-145 


mmu04666:Fc gamma R-mediated 
phagocytosis 


ARF6, CFL2, LAT, ARPC5, CRKL 


5 


0.01 


mmu-miR-145 


mmu04520: Adherens junction 


YES1, ERBB2, CTNND1, ACTG1, SMAD3, ACTB, TGFBR2 


7 


1.18E-06 


mmu-miR-145 


mmu04310: Wnt signaling pathway 


PPP3CA, FZD9, MYC, CTNNBIP1, SMAD3, SENP2, WNT5B 


7 


3.62E-04 
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Table 1 Predicted targets of modulated miRNAs upon RABV infection involved in immune response pathways 

(Continued) 



mmu-miR-145 


mmu04360: Axon guidance 


PPP3CA, CFL2, SRGAP1, DPYSL2, SEMA3, EFNB3, SRGAP2, ABLIM2, 
SEMA6A, PLXNA2, SEMA3D 


11 


0.001 


mmu-miR-145 


mmu041 10: Cell cycle 


MCM5, MYC, SMAD3, CDK6, SFN, 0RC2L 


6 


0.007 


mmu-miR-145 


mmu04350: TGF-beta signaling pathway 


MYC, SMAD3, INHBB, TGFBR2, SMAD5 


5 


6.96E-06 


mmu-miR-145 


mmu04510: Focal adhesion 


ERBB2, ACTG1, FLNB, ACTB 


4 


0.018 



studies have also used profiling technology to evaluate 
the modulations in miRNA expression that occur in re- 
sponse to viral infection. For example, Li and colleagues 
performed miRNA profiling in the lungs of mice 
infected with influenza and found that cellular miRNA 
might be a contributing factor to the extreme virulence 
of the influenza virus [15]. Bhomia and colleagues sug- 
gested that host miRNAs were significantly modulated 
in mouse brains upon VEEV infection [16]. 

Although many computational approaches have been 
developed to predict miRNA targets using sequence in- 
formation, their accuracy is limited [27-30]. To increase 
reliability, we used three web-based target prediction 
databases, TargetScan, MicroCosm, and Targets. The 
pathways listed in Table 1 and Additional file 3: Table S2 
are among the most significant, as indicated by predic- 
tions from all three algorithms. It has been shown that 
miRNA-induced down-regulation of target genes pro- 
vides opportunities to develop new approaches to target 
identification and validation using high-throughput ex- 
pression profiling [31]. Gene expression profiling data 
have been used to identify functional targets of miRNAs 
[32-34]. MiRNAs and the mRNAs that they target for 
degradation can be expected to exhibit an inverse ex- 
pression relationship. Researchers established a strategy 
for miRNA target identification using these inverse rela- 
tionships as predicted from the paired expression profiles 
[35]. In our recently published work, we also identified 
miRNA targets using these inverse relationships as 



predicted from the paired expression profiles [17]. In the 
present study, we simultaneously collected miRNA and 
DNA microarray data from the same samples and com- 
pared the predicted targets of significantly modulated 
miRNAs to the gene expression profiles of RABV-infected 
mouse brains. This showed that some of the predicted 
miRNA targets were correlated with the mRNA expres- 
sion profile. 

Host defense against viral invasion requires induction 
of appropriate innate immune responses. Upon recogni- 
tion of viral components, host cells become activated 
and produce type I IFN and proinflammatory cytokines 
[36,37]. A suitable amount of type I interferon (IFN) 
induces cellular resistance to viral infection and apop- 
tosis of virus-infected cells [38]. However, viruses have 
developed several strategies to evade and subvert the im- 
mune responses mediated by type I IFN, including har- 
nessing host miRNAs. A recent study demonstrated that 
the vesicular stomatitis virus (VSV), family Rhabdoviri- 
dae, can induce up-regulation of miR-146a, which 
feedback-inhibits RIG-I-dependent IFN-I production in 
macrophages [39]. In our study, enrichment of KEGG 
pathways revealed that the predicted target genes of dif- 
ferentially expressed miRNAs upon RABV infection may 
involve Jak-STAT signaling pathway. The Jak-STAT path- 
way is initiated in response to cytokines, such as inter- 
leukins and IFNs, and growth factors. To invasive innate 
immune response of host, RABV interrupts IFN Jak-STAT 
signaling in a manner of activation-dependent targeting of 



Table 2 DNA microarray and qRT-PCR analysis of expression of miRNA targets 



MicroRNA 


Accession 


MiRNA 


qRT-PCR 


Targets of 


Accession 


DNA 


qRT-PCR 




number 


microarray 


(Fold: T or j) 


microRNAs 


Nun 


nber 


microarray 


(Fold: 


Tor |) 






(Fold: T or j) 










(Fold: T or j) 






mmu-miR-691 


MI0004659 


25.29x,| 


16.71 ±3.24x,| 


MAPKAPK3 


NM_ 


J 78907 


1.57x4 


1 0.02 ± 


: 4.1 1x4 


mmu-miR-691 


MI0004659 


25.29x,| 


16.71 ±3.24x,| 


IL25 


NM_ 


_080729 


1.45x # J. 


12.75 ± 


: 2.91x4 


mmu-miR-377 


MI0000794 


12.90x,t 


11.22±2.87x,| 


SOCS4 


NM_ 


_080843 


1.35x4 


9.52 ± 


: 2.56x4 


mmu-miR-1935 


MI0009924 


10.06x,| 


9.94±2.63x,| 


CCR10 


NM_ 


_007721 


1.39x4 


4.58 ± 


: 2.03x4 


mmu-miR-1935 


MI0009924 


10.06x,| 


9.94±2.63x,| 


NFAT5 


NM_ 


J 33957 


1.59x4 


7.69 ± 


: 3.1 7x4 


mmu-miR-203 


MI0000246 


2.79x,| 


3.73 ± 1 .42x,| 


VEGFA 


NM_ 


,001025257 


1.37x4 


4.87 ± 


: 1.96x4 


mmu-miR-290-5p 


MI0000388 


2.44x,| 


6.79 ± 2.3 1x,T 


CD40 


NM_ 


.011611 


1.3x4 


5.43 ± 


: 2.76x4 


mmu-miR-145 


MI0000169 


0.40x,| 


4.82±1.58x,| 


ARF6 


NM_ 


_007481 


1.21x4 


5.38 ± 


: 2.77x4 



Data from qRT-PCR are shown as mean ± standard deviation (SD) of one representative experiment. Similar results were obtained in three independent 
experiments. 



Zhao et al. Virology Journal 201 2, 9:1 59 
http://www.virologyj.eom/content/9/1 /1 59 



Page 9 of 13 



STAT1 and STAT2 [40]. In this study, the target genes of 
modulated miRNAs were found to be involved in Jak- 
STAT signaling, including JAK2, SOCS4, which are targets 
of miR-377; TYK2 and IL12A, targets of miR-691; IL6, a 
target of miR-190; IFNAR2, a target of miR-135a*; 
IFNGR2 and AKT1, targets of miR-290-5p; and SOCS3 
and AKT2, targets of miR-203. This suggests that the Jak- 
STAT pathway may be affected by RABV-inducible cellu- 
lar miRNAs (Table 1). 

Recent studies have revealed the important regulatory 
roles played by cytokines and their receptors in RABV 
infection. One study showed that over-expression of 
cytokine CCL3 (MlP-la) in mouse brains decreased 
RABV pathogenicity [41]. The same research team also 
demonstrated that MlP-la not only reduces viral patho- 
genicity but also enhances immunogenicity by recruiting 
dendritic cells and B cells to the sites of immunization, 
lymph nodes, and blood [42]. We observed that several 
targets of differentially expressed miRNAs are involved in 
cytokine-cytokine receptor interaction. These included 
FASL, IL18RAP, and KTTL, which are targets of miR-377; 
7125, IL12A, TNFSF14, and CLCF1, which are targets of 
miR-691; CCR10, a target of miR-1935; CXCL16, a target 
of miR-190; IL24, IFNG, CXCL16, and CD40LG, targets of 
miR-135a*; IL18R1, CD40, CXCL12, and CSF, targets of 
miR-290-5p; and IL24, XCL1, and IL12B, targets of miR- 
203 (Table 1). Our findings showed that modulated miR- 
NAs may regulate the functions of cytokines during 
RABV infection. 

The MAPK signaling pathway has been shown to regu- 
late the expression of genes involved in the immune re- 
sponse to pathogens. Viral infection can induce activation 
of the MAPK signaling pathway [43,44]. RABV infection 
induces MAPK and NF-kB activation, which have been 
found to regulate chemokine expression in microglial cells 
[45]. Some key MAPK signaling pathway-related target 
genes were identified in the present study. These included 
SPRED2, a target of miR-691, MAP3K12, a target of miR- 
1935, MAPKAPK3 and MAP3K1, targets of miR-290-5p, 
and MAP3K13, MAP4K3, and MAP3K1, targets of miR- 
203 (Table 1). This demonstrated that RABV-induced cel- 
lular miRNAs might be involved in the MAPK signaling 
pathway after RABV infection. 

In summary, the results of the present study provide 
evidence that specific miRNAs are modulated in the 
street-RABV-infected brain. This result was found to be 
completely different from the expression profiles of host 
miRNAs in the CNS of mice infected with the 
laboratory-fixed strains of RABV. Considering that this 
was verified by both DNA microarray and qRT-PCR, we 
suggest that the modulated miRNAs might affect the 
biological processes of cells during RABV infection. Our 
study suggests that host miRNAs might be an important 
class of targets and may play a key role in regulating 



gene expression in response to highly pathogenic RABV 
infection of the CNS. 

Conclusion 

In summary, our findings suggested that street RABV in- 
fection resulted in significant changes in the expression 
of multiple miRNAs in mouse brains. The modulated 
miRNAs might regulate biological processes of cells dur- 
ing RABV infection. The predicted target genes of these 
differentially expressed miRNAs are involved in immune 
responses in the host. MiRNA and mRNA profiles 
obtained in this study might help elucidate the regula- 
tory mechanisms that mediate the host response to 
RABV exposure. 

Methods 

Viruses 

RABV street rabies virus Fujian strain, isolated from a 
rabid dog in Fujian Province, was used for this study. 
Viral stocks were prepared as described elsewhere with 
minor modifications [46]. Briefly, Three-day-old suckling 
mice were intracerebrally (i.e.) infected with 30 ul of 
viral sample. When moribund, mice were euthanized 
and brains were removed. A 10% (wt/vol) suspension 
was prepared by homogenizing the brain in Dulbeccos 
modified Eagles medium (DMEM, Gibco CA, U.S.). The 
homogenate was centrifuged to remove debris, and the 
supernatant was collected and stored at -80°C. The viral 
titers were determined in triplicate on monolayer cul- 
tures of mouse neuroblastoma cell (NA) as described 
previously [47]. 

Animal infection and assessment of clinical signs 

Six-to-eight-week-old female BALB/c mice were obtained 
from the Changchun Institute of Biological Products, 
China. Animals had access to food and water ad libitum. 
All the experiments with live virus challenge were carried 
out at the bio-safety level 2 (BSL-2) facilities of the Key La- 
boratory of Jilin Province for Zoonosis Prevention and 
Control, Institute of Military Veterinary, Academy of Mili- 
tary Medical Sciences. 

The animal experiments were conducted with prior 
approval from the Animal Welfare and Ethics Committee 
of Institute of Military Veterinary, Academy of Military 
Medical Sciences under the permit number (SCXK- 
2002-018). All manipulation of the mice satisfied the 
requirements of the Regulations of Experimental Animal 
Administration of China. 

Mice in the experimental mice were infected with 10 5 
focus-forming units (FFU) of RABV in 30 ul of DMEM 
via the i.e. route. Mice in the control group were mock- 
infected i.e. with DMEM containing uninfected brain 
homogenates for use as controls. Infected animals were 
observed twice daily for 21 days for the development of 
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rabies. Disease progression and mortality were moni- 
tored. Clinical signs were scored as described elsewhere 
using a scale of 0 to 5: 0, no clinical signs; 1, disordered 
movement; 2, ruffled fur, hunched back; 3, trembling 
and shaking; 4, complete loss of motion (complete par- 
alysis); 5, death [48]. 

Tissue collection and total RNA isolation 

Tissue collection and total RNA isolation were performed 
as described elsewhere [18]. Briefly, mice at 7 dpi were 
anesthetized with ketamine-xylazine (1.98 and 0.198 mg 
per mouse, respectively) and euthanized. Brains were har- 
vested and stored in RNAlater (Ambion TX, U.S.) at -80°C 
for total RNA extraction. Total RNAs were isolated from 
the entire homogenized brains using the mirVana™ 
miRNA Isolation Kit (Ambion TX, U.S.). The integrity of 
total RNA was analyzed by Agilent 2100 Bioanalyzer 
(Agilent Technologies, CA, U.S.). 

Taqman PCR quantification of viral loads 

To determine viral load in infected brain tissues, Taq- 
Man qRT-PCR was performed on RNA samples using 
Custom TaqMan® Gene Expression Assays (Applied 
Biosystems, CA, U.S.). The primers specific to the N gene 
of the RABV Fujian strain (forward: 5 -GTGGGTACTGT 
TGTCACTGCTTA-3 ' and reverse: 5 -GTGAGATTTAT 
CTGCTTTATGAACCCTGTA-3 ') and probe (FAM- 
TCCTGAGCAATCTTC-NFQ) and the protocols for 
TaqMan qRT-PCR were used as described by our previous 
study [18]. The TaqMan PCR was performed using Bril- 
liant II qPCR Master Mix (Agilent Technologies, CA, U.S.) 
in an Mx3005P apparatus (Agilent Technologies, CA, U.S.) 
according to the manufacturer s instructions. TaqMan runs 
of experimental samples contained at least three replicates 
with no-template or no-primer controls. Real-time PCR 
was performed in reaction mixtures including 12.5 \A of 2x 
QPCR master mix (Agilent Technologies, CA, U.S.), 1 \A 
of 20x Custom TaqMan® Gene Expression Assays (Ap- 
plied Biosystems, CA, U.S.), 0.375 \A of diluted reference 
dye (Agilent Technologies, CA, U.S.), and nuclease-free 
PCR-grade water to a final volume of 25 ul. The PCR 
conditions were (i) 95°C for 2 minutes and (ii) 40 cycles 
of 95°C for 5 seconds and 60°C for 20 seconds. A standard 
curve was generated from serially diluted RABV N RNAs 
of known copy numbers, and the copy numbers of sam- 
ples were normalized to 1 \ig of total RNA. An absolute 
standard curve method was to calculate the copy numbers 
of RABV N mRNA in mouse brain tissue [49] . To exclude 
contamination of genomic DNA, control cDNA reactions 
in which reverse transcriptase was omitted were prepared 
in parallel as described elsewhere [50]. These were uni- 
formly negative. 



uParaflo miRNA microarray assays 

Three RABV-infected and mock-infected mice at 7 dpi 
were randomly selected for miRNA microarray analysis. 
(iParaflo miRNA microarray assays were outsourced to 
LC Sciences (Houston, TX, U.S.). The assay was per- 
formed on a 5 (ig total RNA sample, which was size- 
fractionated using a YM-100 Microcon centrifugal filter 
(Millipore, MA, U.S.). The small RNAs (<300 nt) isolated 
were 3 '-extended with poly-(A) tails using poly-(A) 
polymerases. An oligonucleotide tag was then ligated to 
each poly- (A) tail for later fluorescent dye staining. 
Hybridization was performed overnight on a [iParaflo 
microfluidic chip using a micro-circulation pump (Atactic 
Technologies, TX, U.S.) [51,52]. The microfluidic chips 
each contained a detection probe consisting of a chem- 
ically modified nucleotide coding segment complementary 
to target microRNA (from miRBase, http://microrna.san- 
ger.ac.uk/sequences/) or other RNA (control or customer 
defined sequences) and a spacer segment of polyethylene 
glycol to extend the coding segment away from the sub- 
strate. The detection probes were made by in situ syn- 
thesis using PGR (photogenerated reagent) chemistry. 
The hybridization melting temperatures were balanced 
by chemical modifications of the detection probes. 
Hybridization used 100 \i\ 6xSSPE buffer (0.90 M NaCl, 
60 mM Na 2 HP0 4 , 6 mM EDTA, pH 6.8) containing 
25% formamide at 34°C. Post-hybridization detection 
used fluorescence labeling with tag-specific Cy5 dyes. 
Hybridization images were collected using a laser scan- 
ner GenePix 4000B (Molecular Devices) and digitized 
using Array-Pro image analysis software (Media Cyber- 
netics, MD, U.S.). Raw data were obtained for further 
analysis. 

DNA microarray assays 

Three mice at 7 dpi were randomly selected from the 
RABV- and mock-infected groups for DNA microarray 
analysis. mRNA microarray assays were outsourced to Phal- 
anx Biotech Group. Inc. (Hsinchu, Taiwan). Fluorescence- 
labeled cRNA was prepared from 5 \xg of total RNA using 
a Message AMPTM aRNA Kit (Ambion, TX, U.S.) and 
Cy5 dye (Amersham Pharmacia, NJ, U.S.). Fluorescent tar- 
gets were hybridized to the Mouse OneArray™ Whole 
Genome DNA microarray (Phalanx, Hsinchu, Taiwan) 
containing 31,802 oligonucleotides, including 29,922 
mouse genome probes, and 1,880 experimental control 
probes. After an overnight hybridization at 50°C, non- 
specific binding targets were washed in three different 
washing steps, and the slides were dried by centrifuga- 
tion and scanned using an GenePix 4000B (Molecular 
Devices, CA, U.S.). The Cy5 fluorescent intensity of each 
spot was analyzed using GenePix 4.1 (Molecular Devices, 
CA, U.S.). Raw data were obtained for further analysis. 



Zhao et al. Virology Journal 201 2, 9:1 59 
http://www.virologyj.eom/content/9/1 /1 59 



Page 11 of 13 



Table 3 Primers for selected genes analyzed using quantitative real-time PCR 

Gene symbol Gene title Forward primer 5'— > 3' Reverse primer 5' ->3' Amplicon (bp) 



MAPKAPK3 


Mitogen-activated protein kinase-activated 
protein kinase 3 


TATOTGTGGCTCCTGAGGTCCT 


TCATOTCCCACCOTAGGT 


112 


IL25 


Interleukin 25 


TGOTGGAGCGCAGGCTCTA 


AGAGAGGGTOGCCCGTA 


113 


S0CS4 


Suppressor of cytokine signaling 4 


ATGTOAAATOCTCTAAGAAG 


CGTAGAACGTCAAAGACAAA^ 


124 


CCR10 


Chemokine (C-C motif) receptor 10 


ACCCAGTGTCTCCCTGATG 


TOCAGTCGGTCCCGGTOA 


136 


NFAT5 


Nuclear factor of activated T-cells 5 


CACCTTCTTCCCCCATTTTAT 


GTCmCAAAAAGGGTOA 


122 


VEGFA 


Vascular endothelial growth factor A 


CCTATOCCCTCTOAATCGT 


AGGAACTGAAGAGAGACCT 


125 


CD40 


CD40 antigen 


TCGGOTOTCTCCAATCAGT 


TGAAAmGGTGTCTACTGT 


126 


ARF6 


ADP-ribosylation factor 6 


GGAAGGTGCTATCCAAGATOT 


TAACCGGGTCGAAGTOAACAT 


108 


GAPDH 


Glyceraldehyde-3-phosphate dehydrogenase 


CTCAACTACATGGTCTACATGTO 


AmGATGTOGTGGGGTCTCGCTC 


142 



Quantitative real-time PCR 

The protocols for qRT-PCR of miRNAs and mRNAs 
were used as described in our previous study [17]. The 
expression of specific miRNAs was analyzed with 
quantitative real-time PCR using a LNA™-based qRT- 
PCR kit, miRCURY LNA™ Universal RT microRNA 
PCR System (Exiqon, Denmark) according to manufac- 
turers instructions. Total RNA from brain tissues of 
RABV-infected and mock-infected mice was converted 
to cDNA using universal reverse transcriptase primers 
(Exiqon, Denmark). cDNA samples were diluted 1:80 
in nuclease-free water and then PCR was amplified 
using SYBR Green Master Mix and specific LNA™ 
miRNA primers (Exiqon, Denmark) for mmu-miR-691 
(target sequence AUUCCUGAAGAGAGGCAGAAAA), 
mmu-miR-377 (target sequence AUC AC AC A A AG G C A 
ACUUUUGU), mmu-miR-1935 (target sequence AGG 
CAGAGGCUGGCGGAUCUCU), mmu-miR-190 (target 
sequence UGAUAUGUUUGAUAUAUUAGGU), mmu- 
miR-203 (target sequence GUGAAAUGUUUAGGACC 
ACUAG), and mmu-miR-145 (target sequence GUCC 
AGUUUUCCCAGGAAUCCCU). Samples were run in 
duplicate on an Mx3005P QPCR system (Agilent Tech- 
nologies, CA, U.S.). The data were normalized to a U6 
RNA control (Exiqon, Denmark) and relative expression 
was calculated using the 2~ AACt method [53]. 

For mRNA quantitative real-time PCR, reverse tran- 
scription was performed using a Reverse Transcription 
System (Promega, WI, U.S.) according to manufacturer s 
instructions. The names of the genes and their primers 
are given in Table 3. Real-time PCR was performed 
using Brilliant® II SYBR® Green QPCR Master Mix 
(Agilent Technologies, CA, U.S.) in an Mx3005P QPCR 
system (Agilent Technologies, CA, U.S.) according to 
the manufacturer s instructions. Expression of the gene of 
interest was normalized to glyceraldehyde-3-phosphate 
dehydrogenase (GAPDH). The expression levels of genes 
were measured in terms of threshold cycle value (CT) 
using the 2" AACt method [53]. 



Microarray data analysis 

Data from microarrays were analyzed as described in 
our previous study [17]. First, raw data from both 
miRNA and DNA microarrays, with replicates, were 
averaged. Probes with median detection P value > 0.05 
and flags reported as absent were filtered out. Expression 
data were normalized using a variance stabilizing 
method VSN in a Bioconductor (http://www.bioconduc- 
tor.org) [54]. Differential expression of miRNAs were 
assessed by eBayes (Empirical Bayes Statistics) from 
LIMMA [55]. Differentially expressed miRNAs with 
P<0.01 were selected for further analysis. Microarray 
data of miRNA and target genes were found to be 
MIAME compliant and have been submitted to the 
Gene Expression Omnibus (GEO) database (http://www. 
ncbi.nlm.nih.gov/geo/; accession number: GSE26269 
(miRNA data); GSE26270 (mRNA data)). 

Target prediction and functional enrichment of 
differentially expressed miRNAs 

Target prediction and functional enrichment of differen- 
tially expressed miRNAs were completed as described in 
our previous study [17]. Three continuously updated 
miRNA target prediction databases, TargetScan, Micro- 
Cosm Targets, and miRanda, were used to infer the targets 
of differentially expressed miRNAs. The intersection of 
these datasets was taken as reliable. Gene ontology (GO) 
category and Kyoto Encyclopedia of Genes and Genomes 
(KEGG) pathway analyses of target genes of differentially 
expressed miRNAs were performed using the web-based 
tool, Database for Annotation, Visualization, and Inte- 
grated Discovery (DAVID, http://david.abcc.ncifcrf.gov/) 19 . 

Additional files 



Additional file 1: Figure SI. Target prediction of differentially 
expressed miRNAs upon RABV infection. Three databases TargetScan, 
MicroCosm Targets, and miRanda, were used to infer the targets of 
differentially expressed miRNAs. 
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Additional file 2: Table SI. Significantly enriched GO terms among 
targets of differentially expressed miRNAs. 

Additional file 3: Table S2. Significantly enriched KEGG pathways of 
differential expressed miRNAs in mouse brain upon RABV infection. 
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